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Abstract 

A binary quenched-annealed hard core mixture is considered in one dimension in order to model 
fluid adsorbates in narrow channels filled with a random matrix. Two different density functional 
approaches are employed to calculate adsorbate bulk properties and interface structure at matrix 
surfaces. The first approach uses Percus' functional for the annealed component and an explicit 
averaging over matrix configurations; this provides numerically exact results for the bulk partition 
coefficient and for inhomogeneous density profiles. The second approach is based on a quenched- 
annealed density functional whose results we find to approximate very well those of the former 
over the full range of possible densities. Furthermore we give a derivation of the underlying replica 
density functional theory. 
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I. INTRODUCTION 

The one- dimensional hard rod model [l| continuous to be an invaluable test bed for theo- 
retical work as it provides the possibility to compare approximations to exact results. Recent 
examples of this strategy include investigations of depletion interactions in binary mixtures 
where one of the components is viewed as an agent that mediates an effective interaction be- 
tween particles of the other component J3 y| ) a model colloid-polymer mixture 4] where par- 
ticles representing polymers can freely penetrate, dynamical density functional theory ^ t \ 
concerned with time-dependent transport phenomena, as well as a model porous medium 7 1 
of lines of random length accessible to the fluid particles. As concerns equilibrium statistical 
mechanics, Percus' exact free energy functional for pure systems J8| and (additive) mix- 
tures of particles with different sizes J3] , provides a framework to compute thermodynamics, 
density distributions, and correlation functions in arbitrary inhomogeneous situations. 

Fluids adsorbed in disordered matrices are often described in the context of so-called 
quenched-annealed (QA) fluid mixtures, where particles of the quenched component act as 
randomly distributed obstacles exerting an external potential on particles of the annealed 
component 10|,lll|- The matrix particles are distributed according to the Hamiltonian of 
the quenched component, and hence can be treated with liquid state theory. The crucial 
difference to an equilibrium system of two annealed components is that the distribution of 
matrix particles is unaffected by the presence of the annealed component. As one is interested 
in the typical behavior of the system, a double average over the annealed degrees of freedom 
and over the quenched disorder is required. A short overview of this theoretical framework 
is given below. One standard approach to tackle fluid structure and phase behavior of QA 
models is via the replica Ornstein-Zernike relations |lOl . Ill| supplemented with appropriate 
closure relations. Many standard liquid integral equation theories have been carried over 
to the QA case. The basic quantities in terms of which these theories are formulated are 
two-body (and possibly higher) correlation functions. 

Recently it was proposed to rather work directly on the level of the free energy functional, 
pre-averaged over the disorder [l2]. Correlation functions can be obtained subsequently, 
in particular the hierarchy of direct correlation functions is obtained through functional 
differentiation with respect to the density fields. The advantage of formulating the theory 
on the one-body level of the density fields is that inhomogeneous situations, e.g. at free 



interfaces or caused by external fields like e.g. gravity, are straightforward to treat. Hence 



all advantages of equilibrium DFT 



m 



IJ] apply to this QA or replica DFT. Also the 
disadvantage applies; in general the density functional is unknown. Moreover, whether 
one can learn anything about the important out-of-equilibrium behavior of QA systems, 
like hysteresis in sorption isotherms, is questionable. Based on Rosenfeld's fundamental- 
measure theory (FMT) for hard sphere mixtures [l5[, and the subsequent discovery that 
one can construct density functionals by imposing the correct behavior upon dimensional 
reduction (i.e. in situations of extreme confinement in one or more spatial directions through 
suitably chosen external potentials) [la, ll2| , this DFT treats QA mixtures with either hard 
or ideal interactions. 

Previous tests of QA DFT include the comparison with results from Monte Carlo (MC) 



simulations for the partial pair correlation functions in hard sphere systems [12| and for den- 
sity profiles across the surface of a porous medium, modeled as a step-like density distribution 
of (freely overlapping) matrix spheres [l^ . Such an interface was also investigated in a model 
of a random fiber network, being represented by quenched configurations of infinitely thin 
needle-like particles. Again comparison with simulation results shows satisfactory agreement 
with DFT results [l9|. Much work has been devoted to a model colloid-polymer mixture 
where the colloids are represented as hard spheres and the polymers as freely overlapping 
spheres. The crucial step beyond hard sphere systems is the occurrence of a fluid-fluid phase 
transition, and the questions how capillary condensation occurs inside a porous medium |20| 
and what the structure of the interface between demixed fluid phases is like J2]| were treated. 
A further, very promising, line of research is the application of the approach to lattice mod- 
els. Note that using lattice models insight into hysteresis behavior and the relation to the 



m 



Combining the QA-DFT 



appearance of a complex free energy landscape was gained 

approach with the very powerful lattice DFT by Lafuente and Cuesta [2j, |25[ , freezing in a 

two-dimensional lattice model was investigated J2a|. 

In this work we consider the one-dimensional hard rod model adsorbed in a quenched 
matrix of rods. Two cases of interactions between the quenched particles are considered. In 
the first case the rods interact with a hard core potential, hence we deal with a binary QA 
hard rod mixture. In the second case the matrix particle are ideal (non-interacting) amongst 
each other, but interact with a hard-core potential with particles of the annealed component; 
this model is the QA analog of the model colloid-polymer mixture of Ref. |^, obtained by 



quenching the polymers. We use two different density functional approaches to tackle the 
properties of these models in both homogeneous and inhomogeneous (on average over disor- 
der) situations. In the first approach the matrix is treated explicitly on the level of particle 
coordinates distributed according to the matrix Hamiltonian, and practically generated with 
a Monte Carlo procedure. For each matrix configuration we use Percus' functional Isj] to 
obtain adsorbate properties, and the disorder- average over matrix configuration is carried 
out numerically by brute force generation of many (of the order of 1000) matrix realizations. 
The second approach is the QA DFT working directly on the level of the disorder-averaged 
adsorbate density profile, which is obtained via minimizing the disorder-averaged grand po- 
tential of the adsorbate, using the matrix density profiles as a fixed input. As the average 
over disorder is taken a priori, the corresponding Euler-Lagrange equation yields directly 
the averaged density profile. We compare results from both theories in bulk via calculation 
of the partition coefficient, which is the ratio of adsorbate density inside the matrix and the 
density in a bulk reservoir that is in chemical equilibrium. The results from the explicit 
matrix averaging procedure, which we also check against an independent elementary calcu- 
lation, agree well with those from the QA DFT over the full range of accessible densities. 
Deviations appear at high densities, which we can trace back to an incorrect behavior of the 
QA-DFT near close-packing. As a generic inhomogeneous situation we consider a surface 
of the model porous matrix, that is generated by a hard wall acting on the matrix particles 
(before the quench) . The wall is then removed and leaves a halfspace of bulk (free of ma- 
trix particles). The adsorbate is found to exhibit density oscillations on both sides of the 
interface, with significantly smaller amplitude than at a hard wall. 

The paper is organized as follows. In Sec. |n] we define the model and give an overview 
of its statistical mechanics and the replica trick (which can be safely skipped by an expert 
reader). Sec. IIIII is devoted to both density functional methods. In Sec. IIVI results are 
presented and we conclude in Sec. El 



II. THE MODEL 

A. Definition of the interactions 

We consider a quenched-annealed fluid mixture of a quencfied species witli Nq particles 
with one-dimensional (Id) position coordinates xi,X2, ■ ■ -Xno and an annealed species 1 
with A^^i particles with Id position coordinates Xi,X2, . . .Xtvi- The particles interact with 
pairwise potentials (for pairs 07 = 00, 01, 11) given by 

00 x< (a« + a^)/2 
(pa^ix) = < (1) 

I otherwise, 

where x is the center-center distance between two particles; cTq is the diameter (length) of 
particles of species a = 0, 1, see Fig. [T^ for an illustration. This describes our first case 
of the hard core matrix. The total potential energy due to particle-particle interactions is 

Voo + K)i + ^11) with contributions 

K)o = ^ 5Z ^ooi\xi-Xj\), (2) 

i=l j=i+l 

No Ni 

^01 = Y.Y.^oi{\x.-X,\), (3) 

i=i j=i 

Ni Ni 

^11 = E E 011(1^^-^.1)- (4) 

We furthermore consider the influence of external potentials, 0^^*(x), acting on species a = 
0, 1, respectively. In particular (j)'^'^{x) acts on particles of species before they are quenched, 
i.e. their density distribution is that generated in response to 0q^*(x). The resulting total 
external potential energy is V^'^^* + ^i^^*, with contributions 

No Ni 

Vo"-' = E C'(^.), vr' = E ^T\x,). (5) 

In our second case we consider ideal (non-interacting) matrix particles, i.e. 

000 (a;) = 0, (6) 

valid for all distances x; the two remaining interactions are unchanged, i.e. 0oi(a;) and 0ii(x) 
are hard core potentials given through (Q); see Fig. ^ for an illustration. The size ratio 



s = o"i/(To is a geometric control parameter. In the numerical results presented below we 
will restrict ourselves to equally sized particles, s = 1, and furthermore to situations where 
0ft (x) = 0. 

B. Partition sum, grand potential and replica trick 

We first make the statistical mechanics of the quenched- annealed mixture explicit. 
For notational convenience the grand canonical trace over matrix coordinates is denoted 
hy J do = Y1'no=o i-^o^-^o°) J dxi . . . J dxN^ and that over adsorbate coordinates by 
J dl = Y^'^ ^g (A'^ilA]^ ^ ) J dXi . . . J dXj\f^ , where A^ is the (irrelevant) thermal wavelength 
of species a = 0,1, and the position integrals run over the total system volume V. The (equi- 
librium) grand partition of the matrix particles under the influence of the external potential 



bn^^(x) is 



o(/io, T,V)= / rfOe-''(^««+^o"-'^°^«), (7) 



where (3 = l/{kBT), ks is the Boltzmann constant, T is absolute temperature and /ij is the 
chemical potential of species i = 0,1. The grand potential for the matrix is then 

Qo{fio,T,V) = -(3''lnEo{fio,T,V). (8) 

For fixed matrix configuration {xj} the grand potential of the adsorbate is 

fii({x,},/ii,T,r) = -/?-Mn /'rfle-^(^"+^«i+^r*-^i^i), (9) 

depending explicitly on {xi} through Vqi, see (jHI). Note that from the viewpoint of the 
1-particles, Vqi + V^^^ is the total external potential energy. A (grand canonical) average 
over matrix configurations yields the disorder-averaged grand potential of the adsorbate, 

fii(/io,/xi,r, V) = So-i(^o,r, V)y"t/0e-^(^"+^o-'-MoA^o)^^(|3,.}^^^^y^V^). (10) 

Note that (fTUj) has a different structure than that of the partition sum of an equilibrium 
mixture due to the appearance of the logarithm, upon inserting Q into (|iup. inside the trace 
over the 0-particles. However, a relation to a multi-component mixture can be established 
using the replica trick: One introduces replicas as s copies of species 1: (paa{x) = 0ii(x), 
for 1 < a < s, where s is an integer. Particles from different replicas are non-interacting, 
(pa-y{x) = for all X and a 7^ 7, but they interact with matrix particles in the same fashion, 



(poa{x) = (f)Qi{x), for 1 < a < s. Then the (equihbrium) partition sum for this (s + 1)- 
component mixture can be written as 

and the grand potential is 

nifxo,fii,T,V;s) = -p-HnE. (12) 

Via analytical continuation in s, and noting that linis^o dx^/ds = \nx, the disorder-averaged 
grand potential, (fTUI). is obtained from the equilibrium grand potential of the replicated 

system as 

fii(/io, /xi, T, V) = lim — fi(/io, /^i, T, V; s), (13) 

s^o ds 

establishing a practical route to tackle the QA system via the replicated equilibrium system. 

III. DENSITY FUNCTIONAL APPROACHES 

The following subsections IIII A[ IIII Bl are valid in arbitrary space dimension d upon trivial 
alterations: spatial integrations become c?-dimensional integrals, hence J dx is to be replaced 
with j d'^x and factors A^ are to be replaced with A^. Although in our subsequent study 
we only deal with hard core interactions, where the dependence on temperature is trivial, 
the formalism also applies to thermal systems. Hence we consider a general binary QA 
mixture with arbitrary pair potentials (pa-ri^), a,7 = 0, 1, not necessarily given through ((H), 
at temperature T inside a volume V. 

A. Equilibrium case 

In the equilibrium DFT formalism, applied to the present case, where both an explicit 
external potential, 0i(x), and the (random) influence of the matrix particles at positions 
{xi} acts on the fluid, the grand potential of the adsorbate component is expressed as a 
functional of its one-body density distribution. 



fii({x,}, [pi],/ii,T, V^) = ^'^[pi] +^<=-[pi 



+ dxpi{x) J0f*(x) + ^0oi(x-a;i) j -/i 



(14) 



where J^"^[pi] = P~^ f dxpi{x)[\n{pi{x)Ai) — 1] is the (Helmholtz) free energy functional of 
the ideal gas, J^'^™[pi] is the excess (over ideal) contribution that arises from interactions 
between (adsorbate) particles, and the term in round brackets is the total external potential 
acting on the adsorbate stemming from the explicit (non-random) external potential, 0^^*(x), 
and the sum over interactions with (randomly distributed) matrix particles. The latter 
contribution is parameterized by the set of matrix coordinates {xj}, and hence Qi depends 
explicitly on {xj}, which we stress in the notation of the l.h.s. of ()14|) . The free energy 
functional JF"^ and JF'^'"^ depend on T and V; this is suppressed in the notation in (fT^ and 
in the following. The minimization condition is 

6n{{xi},[pi],ni,T,V) 



= 0, (15) 

Pi=Pi{{xi},x) 



Spi 

where pi{{xi},x) is the adsorbate density distribution that solves (fTKj) . The value of the 
grand potential is then obtained by reinserting the solution into the grand potential func- 
tional, 

n{{x,}, pu T, V) = n{{x,}, [pi(x, {Xi})],pu T, V), (16) 

from which the average over the disorder, Q{po, pi,T,V), can be obtained via (|Tn|l . In a 

similar way as for the grand potential, the matrix-averaged adsorbate density profile is given 

as 

pi(x) = Ho-i(po,T,\/)y"ciOe-^(^-+^o-'-MoA^o)^^(^^|^.|). (17) 

Note that explicit averages over the matrix configurations are to be performed in (jiup and 
(fT7|) . In the numerical procedure described below, we will carry this out numerically via a 
Monte Carlo procedure. 

B. Quenched- annealed case 

Here we first formulate the equilibrium DFT of the replicated model from which we obtain, 
in the appropriate limit of vanishing number of components, the minimization condition of 
DFT for one quenched and one annealed component. Explicitly assuming absence of replica 
symmetry breaking, hence pi(x) = pa{x),l < a < s, the grand potential functional for the 
replicated equilibrium mixture reduces to 

ni[po,pi],po,fti,T,V;s) = ni[{p^}],{p^},T,V) (18) 



The minimization conditions are 



0, a = 0, 1, (19) 

Opa{x) 

which are two coupled equations for the two unknown functions po{x) and pi{x). At the 
minimum the value of the functional is the true grand potential 

fi(/io, /ii, V, T; s) = Cl{[po, Pi], /io, Pi, T, V; s). (20) 

Via analytic continuation, and Taylor expanding in s around s = 0, one obtains 

n{[po,pi],p.o,l^i,T,V;s) = no{[po],p.o,T,V) + s^idpo, pi], ni,T,V) + 0(s^), (21) 

where ^o is the grand potential of the pure system of 0-particles, that may formally 
be written as Qo{[po], po,T,V) = r2([po, 0],po,Pi ^ —oo,T,V), furthermore Qi[po,pi] = 
lims^o (^^([po. Pi], Po, Pi, 7", V^; s)/t/s. Note that via this definition Qi is also the disorder- 
averaged grand potential, given in (fT^ . 

We decompose both contributions in the standard way: 

fio([po], Po, T, V) = r^lpo] + JTIpo] + J dxpoix) (0^o'*(x) - po) , (22) 

fii([po,Pi],Pi,T,\/) = r''[pi]+J^r{po,Pi] + jdxp,{x){4>l^\x)-fi,), (23) 

which can be viewed as definitions for the excess (over ideal gas) free energy functionals, 
^q^'^IPo] and J-'^^'^[po, pi], that arise from interactions between particles. Note, however, that 
neither jF"^[po] nor a contribution involving 0q^*(x) appear on the r.h.s. of ()23|1 . in contrast 
to the binary equilibrium case. 

We insert the small- s-expansion of the grand potential functional, (j2^, into the equilib- 
rium minimization condition, (fT^ . Performing the limit s — »• 0, the minimization conditions 
for the QA mixture result as 

(^^o([Po],Po,r, V) 



Spo{x) 
5^i([po,pi],pi,T, 1/) 



0, (24) 

0. (25) 



6pi{x) 

The derivation of (|24j) is straightforward. To obtain (j25|) one sets a = 1 in (|19p. i.e. 
differentiates w.r.t. to pi(a;) and divides the resulting equation by s, assuming s > 0, before 
taking the limit s — » 0. Note that (j^ is decoupled from (J2SI), hence in particular po{x) is 
solely determined through ()24|) . The result then serves as an input to ()25|). which is solely 
to be solved for pi(a;). 



C. Excess free energy functionals 

The theoretical approaches described so far incorporate the complexity of the problem i) 
in the excess free energy functionals for the adsorbate, jF'^^^[pi], and an explicit average over 
matrix configurations, and ii) in excess free energy functionals for the matrix component, 
^o"^[Po], and for the adsorbate in the presence of the matrix, jF{'™[po,pi]. For the present 
Id model, the situations is fortunate, as an exact result for J^'^^'^[pi] and JFq'"^[po] is available, 
namely Percus' free energy functional for Id hard rods [8|, |9|. For a one-component system 
of hard rods of species a one writes (using Rosenfeld's terminology) 

^-IpJ = J dxn(^\x)Kc {n'^^Hx)) , (26) 

where the weighted densities, Ua (x) and ria (x), are obtained from the bare density profile 
(of species a) via 



"X + B.a 



n):'{x) = [pa{x - Ra) + Pa{x + Ra)]/2, (27) 



n^^\x) = / dx'p^ix'), (28) 

J x—Bct 

where Ra = o"q/2 is the particle "radius" of species a = 0, 1, and the upper index z/ = 0, 1 
of the weighted density is related to its dimension, which is (length)'""'', where d = 1 is the 
space dimension. The prime in (j26|) denotes differentiation w.r.t. the argument, and $hc(^) 
is the zero-dimensional free energy of hard core particles |la], given by 

^i,,{r]) = {l-r])\n{l-r])+r]. (29) 



Our (approximate) QA functional has very similar structure 12|- We start from the 
generalization of ()26|) to binary mixtures, which is 

•^r[Po,Pi] = [dx J2 n^a\^)^a (4^^(x),nf)(x)) , (30) 

"^ 0=0,1 

where the weighted densities are still given through (jTTj) and (J2HI); and derivatives of the 
zero-dimensional free energy, $, are defined as ^aiVo^Vi) = d^{VQiVi)/9f]a- The particular 
form (j3(jp ensures that the exact result for $ is recovered if the functional is applied to a 
zero-dimensional density distribution, defined as Pi{x) = r]iS{x) for i = 0,1. Hence it can 
be shown via elementary calculation that $(770,771) = J-'i^^[r]o6{x),rii6{x)]. The exact result 
for $(770,^1) is obtained by solving the zero-dimensional limit (where all particles present in 

10 



the system overlap); a detailed calculation can be found in Ref. [l^]. Such a situation can 
be enforced by appropriate external potentials (j)^^{x) = oo if |a;| > L and zero otherwise, 
and represents a "cavity" of size L, where L -C niin(cro, o"i)/2. Hence any two particles 
(of species a and 7) present in the system will overlap, i.e. x < cXa-y, where x is the center- 
center distance between both particles. (As a consequence the density profiles vanish outside 
the cavity, Pa(|x| > L/2) = 0.) This constitutes the crucial simplification that allows to 
calculate the free energy; its excess contribution is independent of L 'l2|. In the case of the 
hard core matrix the result is 

*('7o, ^71) = (1 - r^o - Vi) ln(l - ^0 - Vi) + r^i - (1 - rjo) ln(l - r/o). (31) 

It is interesting to note that in this special case there is a simple relation to the corresponding 
fully annealed binary hard core mixture: $(?70) ^71) = '^hc('7o + ^1) ~ '^hc(^o)- In the case of 
the ideal matrix the result for the excess free energy is 

$(r?o, r^i) = (e-^o - r^i) In (e""" - r],) + rji + rjoe-'"> . (32) 

This does not have a similar relation as above to the corresponding fully annealed binary 
mixture (the colloid-polymer mixture of Ref. |j4|). As a further aside, we note that setting 
'^('705 ^71) = '^hc('7o + ^71) in (jHUj) gives the exact excess free energy functional for equilibrium 
(where both species are annealed) binary hard rods J8|, |9| . 

The density functional thus defined is exact on the second virial level, which can be 
seen by Taylor expanding ^{i]Q,rii) in both arguments to second order, and exploiting the 
property of the weight functions to recover the Mayer bond upon convolution. 

D. Numerical procedure 

In order to calculate density profiles, the minimization is done with a standard iteration 
technique. We chose a very fine grid with spacing 0.0002cr to discretize the space coordi- 
nate. The system is assumed to be periodic in x with length 20a. To generate the matrix 
configurations {xj}, a Monte Carlo scheme is used. We start from an initial configuration 
where the (matrix) particles have equally-spaced positions. Then particle displacements are 
performed according to the Metropolis algorithm, i.e. a new position is only accepted if no 
overlap with any other (matrix) particle or with the wall (introduced below) occurs. 1000 

11 



MC moves per particle are performed for equilibration, and 100 MC moves per particle are 
performed between configurations that are used for data production. For given matrix con- 
figuration {xi} the density profile for the adsorbate component, pi(x, {xj}), is then obtained 
from solution of the minimization condition of the equilibrium DFT, (jl5|) . Results from 
1000 independent matrix realizations for each Nq in the range A^^i = — 20 are then used to 
carry out the average over the disorder and obtain pi{x) via (fTTj). 

IV. RESULTS 
A. Bulk 

The bulk case constitutes the basis of the subsequent interface study, and is considered 
for the following three reasons: i) to assess the accuracy of the QA-DFT, ii) to demonstrate 
the correctness of the matrix-averaging procedure via comparing with an independent el- 
ementary calculation, and iii) to study the (exact) partition coefficient, which we find to 
possess (very unusual) non-monotonic behavior. 

The matrix particles are distributed homogeneously on average, and hence are charac- 
terized by a one-body density distribution po{x) = tiq/ctq = const. As a consequence, the 
adsorbate density distribution is on average pi(x) = r/i/ai = const. We imagine the system 
to be in chemical equilibrium with a reservoir of hard rods of species 1 of packing fraction 
rjl = PiCTi, where pi is the number density in the reservoir; there are no quenched matrix 
particles in the reservoir. The reservoir density sets the chemical potential of 1-particles via 
the (well-known) hard rod equation of state, 

/3/ii = ln(r/[AiM) - ln(l - r^[) + -^. (33) 

l-r]l 

The central quantity that we use to characterize the coupled systems is the partition co- 
efficient K = rji/rjl = Pi/pI, which we will study as a function of the reservoir packing 
fraction rjl. 

Before applying both DFT methods to this case, we first seek^to obtain a benchmark 
result for K{r]l) from an independent, elementary calculation; see 27| for more mathematical 
background and ^] for an alternative application. We start from the probability W{x)dx that 
the nearest-neighbor distance for 0-particles (measured between their centers) is between x 



12 



and X + dx (we call this a gap of size x) , given by 



l-?70 ^ 



^(1-x/ao)) x>ao 



W{x)a^= { '-"^ V^-^'° ^ (34) 

otherwise, 

from which the average nearest-neighbor distance between matrix particles follows as x = 
J dxxW{x) = (Jo/rjo, which is an expected result. To derive (jHH), note that the occurrence 
of a gap of size x is proportional to the Boltzmann weight of the reversible mechanical work 
to create it, hence W{x) oc exp(— Pq^^), where Pq is the pressure of the bulk (matrix) system, 
given for one-dimensional hard rods by ctq/SPo = ^70/(1 " Vo)- ^^^ obtains the precise form 
of ()34p by taking into account the correct normalization, J dxW{x) = 1. 

The average number of 1-particles in a gap of (fixed) size x between two 0-particles is 
iVi(x, /ii) = P~^d\nEi/dfii, where the partition sum of 1-particles in the gap of volume 
X — ao — ai (being accessible to the centers of 1-particles) is 

Si(x,/xi)= [dlexpi-PiVu-^^lN^))= Y. ^-^^^^^^^{x - ao - N,a^f\ (35) 

The mean number of 1-particles in the gap between two neighboring 0-particles, averaged 
over all sizes of the gap, is iV'i(/ii) = J dxW{x)Ni{x, fii). Then the average density of 1- 
particles is pi = Ni{jj,i)/x = Ni{fii)rio/(TQ, from which the partition coefficient results as 
K = pi/p'i = Niai/{xril). Putting things together, 

where the equation of state in the reservoir of 1-particles, (jSBI), can be used to obtain /xi in 
terms of ?][ in ()36|). hence K is solely a function of the packing fractions ?7o, t]1 and of the size 
ratio cri/cTQ. In general, we solve ()36p numerically; we can, however, obtain analytic results 
in both (extreme) cases of high and low reservoir density, 

K{r][ ^ 0) = (1 - Vo) exp(-sr/o/(l - Vo)), (37) 

Kivl - 1) = , ,.7° ,, , , (38) 

exp(sr/o/(l-^o)) -1 

where s = o-i/ctq- ()38|) is obtained by noting that for fii ^ 00 the mean number of 1-particles 
equals the maximal possible number, i.e. Ni{x, /ii -^ 00) = fioor((a; — cro)/o"i), where floor(a;) 
gives the next integer smaller than x. 

13 



In the case of the of the non-interacting matrix, the equation of state of the 0-particles 
is that of an ideal gas, aoPP = rjQ. The gap size distribution is 

W{x) = {r]o/(To)exp{-r]ox/cro), x > 0, (39) 

from which K can be derived following the same steps as above. Again in the two limiting 
cases of either high or low reservoir packing fraction analytic expression can be obtained for 
the partition coefficient, 

K{r][^0) = exp(-(l + s)r/o), (40) 

Kivl - 1) = '-H^^^tnA. (41) 

exp{sr]o) - 1 

In order to obtain K from the QA DFT, we use the predicted equation of state, which is 
in the case of the hard core matrix 

/3/ii = ln(r7iAi/ai) - ln(l - rjo - Vi) + "^^ ^ '"^^ , (42) 

1 - ^70 - ^1 

which we set equal to the chemical potential in the reservoir, (jH^ . For small ri[ this gives 

the exact result, (J37jl . In the high density limit, however, we do not recover (J38|) . but 

obtain K{ril ^ 1) = 1 — rjo- Note that the small r^o-expansion of the exact result, (jHEj), 

is K{rii = 1) = 1 — (1 + s/2)?7o + 0(?7q), differing already in linear oder. For intermediate 

values of r]{ we obtain numerical solutions. 

Results for K a.s a. function of reservoir packing fraction, rjl, are displayed in Fig. |2K for 

three different matrix packing fractions, r^o = 0.1,0.5,0.7. For the lowest value considered, 

t]q = 0.1, K slightly increases as a function of r]'[. Remarkably a maximum is reached at 

ri[ ~ 0.9, significantly smaller than the close-packing limit, ri[ = 1. We could not obtain 

numerical results for 0.9 < 77^ < 1, but the available data clearly tend towards the limiting 

value, obtained through (jHHj) . For the intermediate value r/o = 0.5 the variation of K with 

rjl is stronger and the maximum is more pronounced. For r][ = 0.7 hardly any adsorbate 

particles can enter the matrix and the resulting values are i^ < 0.1, the variation with rjl 

being similar to the above cases. The exact results plotted in Fig. El are obtained from 

the elementary calculation above (lines) and from the DFT with explicit matrix averaging 

described in Sec. IIIIAI Both agree with high numerical accuracy. The result from the QA 

DFT is exact in the low-density limit and stays accurate for intermediate rjl. At r7[ ~ 0.7 

deviations emerge, overestimating the values of K. Moreover, the maximum is absent, 
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and monotonic increase with rjl is found. We hence conclude that the overall performance 

of the QA DFT is very satisfactory, but that the intriguing non-monotonic behavior near 

close-packing is missed. 

The situation in the case of the ideal matrix is similar. The equation of state from QA 

DFT is 

/3/zi = ln(r^iAi/ai) - In (e'^o - r^i) + '"^"'"f + '^\ (43) 

e ^1° — rji 

Again the result from QA DFT for rjl —>■ equals the exact result, (PT|) . In the limit 
?7[ — > 1 the result from QA DFT is K{ril ^ 1) = exp(— ?7o), which overestimates the exact 
expression, (j4ip . See Fig. |21b) for numerical results for i^ as a function of rjl for matrix 
packing fractions rjo = 0.1,0.5, 1. Compared to the case of the hard core matrix, at equal 
densities K is larger, clearly due to the more open void structure of the ideal matrix particles. 
Again K displays a maximum near rjl ~ 0.8, which is not captured within the QA DFT. 
Nevertheless, the overall agreement is again very reasonable, and gives us confidence to turn 
to inhomogeneous situations. 

B. Behavior at the matrix surface 

We consider matrix distributions that are generated by a hard wall, described by the 
external potential 

{0 a; > 
(44) 
oo otherwise, 

acting before the quench. In order to obtain po{x), we solve the minimization condition, 

(J2D), using Percus' functional, given through (j^HJ-iEni), as JFq™ in (j221)- This provides an 

exact (numerical) solution, see Fig. IHt) for matrix density profiles for three different matrix 

packing fractions rjo = 0.1,0.5,0.7. As rjo increases the contact value at the wall increases 

and the layering near the wall becomes more pronounced, i.e. the amplitude of the density 

oscillations grows. 

The density profile of the adsorbate being exposed to such an inhomogeneous matrix 

(note that the hard wall, (jHj), only acts on the matrix particles) are obtained using either 

the explicit matrix-averaging procedure of Sec. IIIIAI or the QA DFT of Sec. IIIIBI In Fig. 

Eb-d we display results for pi(x) for rjQ = 0.1, 0.5, 0.7. For x > far away from the surface, 

the reservoir value is reached, pi{x —^ oo) = p^; for x < deep inside the matrix, the 
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asymptotic value is pi(x -^ — oo) = Kp'l, as studied above. For the lowest adsorbate density 
considered, result for /3/ii = —5 are displayed in Fig. El the crossover between the limiting 
cases happens in a narrow interval —a < x < a. Already for x > a the profile is flat. 
Inside the matrix, however, for x < a there are small oscillations with wavelength of the 
order of a that decay rapidly with increasing distance from the surface. These oscillations 
are not due to the correlations between adsorbate particles, but are merely "imprinted" by 
the inhomogeneous matrix profiles as shown in Fig. ^i. Increasing the adsorbate chemical 
potential, see Fig. Et for /?/ii = 0, leads to stronger structuring outside the matrix, x > a. 
This layering is similar to that at a hard wall, but with significantly smaller amplitude. 
Clearly, this "washing out" is due to the average over the disorder. Note that for a given 
matrix configuration {xi} the matrix particle closest to the surface, say x^^, exerts a hard 
interaction on the fluid particles with Xj > xnq, and indeed the resulting pi(x) is that of 
a hard wall. The disorder-average then leads to the "washing out". For (3fii = 3, shown 
in Fig. IHli, even stronger layering is observed. As expected from the bulk analysis above, 
the QA DFT overestimates K and hence the adsorbate density inside the matrix, which is 
clearly visible for x < —o. The shape of the curves, however, is predicted very accurately. 
Moreover, for x > —a, the QA DFT lies practically on top of the exact result. 

In the case of the ideal matrix particles exposed to the hard wall, Eq. (jUJ, leads to step 
function density profiles, see Fig. |^, hence is very different from highly structured profile 
of the hard core matrix. Fig. IHl The general trends upon varying the matrix density and 
the adsorbate chemical potential, see Fig. Hb-d for results for the same values of /?/ii as 
above, is similar. The ideal matrix allows to consider higher matrix packing fractions, and 
we have gone up to ?7o = 1. It is to be noted that oscillations do appear for x < a, although 
the matrix profile is uniform. These oscillations clearly arise from packing effects between 
adsorbate particles. Again the results from the QA DFT agree very well with those of the 
exact calculation. 

V. CONCLUSIONS 

In conclusion we have applied a recent DFT for adsorbate fluids in random matrices 
to the one-dimensional hard core model. Comparing with analytic and numerical exact 
solutions in bulk and at matrix interfaces we have tested the accuracy of the QA DFT. 
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The general performance is remarkable, both for the bulk partition coefficient and for the 
inhomogeneous density profiles. Subtleties like non-monotonic variation of the partition 
coefficient upon increasing adsorbate density are not captured by the QA DFT. 

We note that by confining mesoscopic colloidal particles in narrow channels J28| one- 
dimensional model fiuids are experimentally accessible. In principle we could imagine mod- 
ifying setups of confining groves as described in Ref. |2^ in order to prepare model systems 
that resemble the model described in the current work. Possible future work could be devoted 
to the impact of quenched disorder on three-dimensional narrow channels 29|, |30| . Further- 
more, it would be interesting to study the asymptotic decay of correlation functions jSll |22| 
of fluids with quenched disorder more systematically. Finally, whether non-monotonic vari- 
ation of the partition coefficient appears in 3d hard core models is an intriguing question. 
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FIG. 1: Model of annealed hard rods in a matrix of quenched rods in one dimension. The quenched 
rods interact with a) a hard core potential, and b) are ideal. 
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FIG. 2: Partition coefficient K = rji/r]\ for Id hard rods immersed in a Id matrix of quenched 
rods as a function of the packing fraction in a reservoir of rods, ij^. Shown are results from the 
elementary calculation, H36|) . (lines and full squares), DFT with explicit matrix averaging (crosses) 
and QA DFT (open squares), a) Matrix particles are interacting with a hard core potential and 
possess packing fraction r/o = 0.1,0.5,0.7 (from top to bottom), b) Matrix particles are ideal and 
posses packing fraction % = 0.1,0.5, 1 (from top to bottom). 
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FIG. 3: Behavior of adsorbate particles near the surface of a matrix of quenched rods immersed in 
a hard core matrix, a) Density profile of the matrix particles, apQ{x), as a function of the scaled 
distance x/a for packing fractions r/o = 0.1,0.5,0.7 (from bottom to top). Adsorbate density 
profiles api{x) are shown for r]Q (from top to bottom) for Pui = — 5 (b), (c), and 3 (d). Results 
from the QA DFT (dashed lines) are compared with those of the exact treatment (full lines). 
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FIG. 4: Same as Fig. |21 but for ideal matrix particles of packing fraction tjq = 0.1,0.5, 1. 
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